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Inelastic transport of electrons through a two- impurity chain is studied theoretically with account 
of intersite Coulomb interaction, U. Both limits of ohmic transport (at low bias) and strongly non- 
ohmic transport (at high bias) are considered. We demonstrate that correlations, induced by a 
finite U, in conjunction with conventional Hubbard correlations, give rise to a distinct transport 
regime, with current governed by two-electron hops. This regime realizes when a single-electron hop 
onto the chain and a single-electron hop out of the chain are both "blocked" due to the finite U, 
so that conventional correlated sequential transport is impossible. The regime of two-electron hops 
manifests itself in the form of an additional step in the current- voltage characteristics, I{V). 

PACS numbers: 72.20.Ee, 71.23.An, 71.55.Jv 



I. INTRODUCTION 

A model of electron transport via two localized states 
positioned in sequence has been studied theoretically for 
more than two decades 1-8 . There are several reasons 
why this model attracts the attention of researchers: (i) 
the model is tractable analytically; (ii) it yields nontriv- 
ial predictions that are in agreement with experimental 
observations 9-15 ; (Hi) it captures physics that is much 
richer than tunnel transport through a single site. In 
particular, the two-site model allows one to reveal vari- 
ous aspects of correlated transport, that are relevant for 
the bulk systems, and to treat correlated transport ex- 
actly. 

In general, for the bulk hopping conductivity, two 
mechanisms are known to cause correlation in the time 
moments of different hops, thus giving rise to the corre- 
lated character of transport. 

The first mechanism is the on-site Hubbard repulsion, 
Uq. In the limit Uq — > 00, the occupancy of a given site 
is restricted to the values of or 1 . The resulting corre- 
lation (dubbed as Hubbard correlation 16 ) reduces to the 
requirement that electron hops occur only between singly 
occupied and empty sites. 

Calculation of the current through a two-site chain 
with Hubbard correlations was carried out in Refs. 1, 4, 
5, 7 both for ohmic 1 ' 4 and strongly non-ohmic 5,7 regimes. 
These calculations, being exact, allow one to test the ap- 
plicability of a standard mean-field description 16 of the 
hopping transport. Within this description, the time mo- 
ments of different hops are assumed completely uncorre- 
cted, while the probability of a given hop is determined 
by average occupations of the constituting sites. 

For finite Uq, double occupation of a single site be- 
comes possible. Together with suppression of the Hub- 
bard correlations, finite Uq makes the transport through 
a chain spin-dependent 6,11 . This is because the second 
electron hopping on a given singly occupied site must, by 
virtue of the Pauli principle, have the opposite spin. 

The second mechanism, that causes correlation be- 



tween different hops is the intersite Coulomb interaction, 
U . Due two a finite U, the energy position of a given site, 
and, thus, the probability of an electron hop onto this 
site, depends on the occupation of the neighboring sites. 
The energy U can be viewed as an additional charging 
energy required for an electron to hop onto and out of 
the two-site chain. 

On average, such a charging, being analogous to the 
Coulomb blockade, impedes the hopping transport. For 
bulk systems this mechanism is commonly accounted 
for within the mean-field theory 16 . The key assump- 
tion adopted to incorporate charging into the standard 
scheme of calculation of the hopping conductivity is that 
this charging amounts exclusively to the depletion of the 
density of states near the Fermi level 17 . Under this as- 
sumption, strongly correlated time evolution of the popu- 
lations of all sites is replaced by completely uncorrelated 
evolution of populations of much fewer sites, the argu- 
ment used for such a replacement being that only the 
hops between these fewer sites govern the transport. 

Another effect of the intersite Coulomb interaction, 
which had received much less attention, is that this inter- 
action opens the possibility for many- electron hops 18 . In 
the course of such hops, two (or more) electrons change 
their spatial positions upon absorption (or emission) of 
a single phonon. This process is analogous to the light 
absorption in the helium atom 19 , in course of which two 
electrons can be excited by one photon. Obviously, such 
an absorption is possible only due to the electron-electron 
interactions. It is also apparent that, as a new transport 
channel, many-electron hops facilitate the transport. 

In the early paper Ref. 18, where many-electron hops 
were first treated analytically, the roles of the two elec- 
trons, participating in the two-electron transition, were 
very dissimilar. While the transfer of the first elec- 
tron occurred between the sites belonging to the current- 
carrying network, the initial and final states of the sec- 
ond electron did not belong to the network. Conse- 
quently, facilitation of transport by two-electron tran- 
sitions amounted to the effective reduction of the activa- 
tion energies of single-electron hops within the current - 
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carrying network. 

Later, on the basis of numerical simulations 20,21 it was 
concluded that many-electron hops may constitute a sig- 
nificant portion of the current-carrying path through the 
sample. However, no analytical theory of hopping trans- 
port with account of many-electron transitions has been 
developed so far. This is because all existing theories are 
based on introducing effective resistors between the pairs 
of sites, which is impossible in the presence of many- 
electron transitions. 

In the present paper, we take advantage of the fact 
that the two-site model is exactly solvable and thus al- 
lows one to study the interplay of all three correlation ef- 
fects, namely, Hubbard correlations, Coulomb-blockade- 
induced correlations, and correlated two-electron hops. 
The fact that the first two correlation mechanisms im- 
pede the transport, while the third one facilitates it, sug- 
gests that such an interplay is nontrivial. 

Our main result is a demonstration of a distinct regime 
of transport through the two-site chain, in which two- 
electron hops dominate the passage of current. This 
regime becomes possible due to intersite interaction, 
U, when both manifestations of this interaction are 
at work: (i) finite U allows two-electron transitions, 
and (ii) it blocks both individual single-electron tran- 
sitions onto and out of the chain. Conceptually, this 
interaction-dominated transport regime is analogous to 
inelastic cotunncling through a Coulomb blockaded quan- 
tum dot 22-24 . However, in terms of its manifestations, 
the important difference between the dot and the two- 
site chain is that the dot contains many levels, so that 
the cotunneling rate is a sum of many contributions. By 
contrast, the transport through a blocked two-site chain 
is governed by a single two-electron hop, leading to dis- 
tinctive temperature and bias dependencies of current. 

We first identify the regime of two-electron hops for 
strongly non-ohmic transport at zero temperature, where 
it manifests itself in the form of additional steps in the 
current- voltage characteristics, I(V). Then we demon- 
strate the relevance of two-electron hops for the ohmic 
transport. 



Similarly, the state, p, in the right lead is coupled to site 
2 with amplitude, ti,p- Then the waiting time, t x , for an 
electron in the left lead to tunnel onto site 1 is given by 
r^ 1 = fi(ei)Ti(si)/h, where the tunneling width, ^(e), 
is defined as 

Ti(e) = 2?r^ |ii ifc | 2 <5(e - £ k ) oc exp 



D - 2di 



(1) 



The corresponding expression for the waiting time, r 3 , 
for electron on site 2 to escape into the right lead reads 



l-/i(e 2 )]r r (e 2 )/ft, with 



r r (e) = 27r V \t2, P \ 2 S(s - £ P ) oc exp 



D~2d 2 



(2) 



Here a is the localization radius of the on-site wave func- 
tions; fi, f r are the Fermi distribution functions in the 
left and right leads, respectively, which determine the 
temperature dependencies of t x and r 3 . 

Transition 1 — > 2 involves tunneling accompanied by 
the emission of a phonon with energy E\ — s%. Thus, 
the temperature dependence of the corresponding time, 
t 2 cx exp[2(c?i + G?2)/a], is weak. 
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II. HOPPING TIMES 

The two-site model is illustrated in Fig. 1. Under the 
applied bias, V , the Fermi levels in the left and right 
leads are shifted V/2 and —V/2, respectively. The two 
sites, 1 and 2, are located at distances d\ and d 2 with 
respect to the center of the barrier of thickness, D. We 
adopt the definition of the energy position of site 1 to be 
ei, when site 2 is empty. Consequently, when site 2 is 
occupied, the energy position of site 1 is £i + U. Anal- 
ogously, £2 and £2 + U are the energy positions of site 2 
for empty and occupied site 1, respectively. We assume 
that there is tunnel coupling with amplitude, £x,k) be- 
tween site 1 and the extended state, k, in the left lead. 



FIG. 1. Schematic illustration of the two-site model. Evo- 
lution of the occupation numbers of sites in the sequential 
hopping regime for blocked (V < V c ) and unblocked (V > V c ) 
transition I — * 1. 

In all previous analytical treatments 1-8 of the two-site 
model it was assumed that the passage of current is gov- 
erned by only three transitions: I — ► 1, 1 — ► 2, and 2 — » r. 
Our key point in this study is that the tunnel couplings 
tj. it and t%„ together with finite intersite interaction, U, 
allow for the additional transition, (1,2) — » (l,r), that 
involves simultaneous change of occupation of both sites. 
Moreover, we will demonstrate that, within a certain bias 
range, these two-electron transitions dominate the pas- 
sage of current through the two-site chain. Calculation 
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of the waiting time, r c , for the two-electron transition is 
similar to the calculation 2 of the inelastic cotunneling 
rate through a quantum dot 
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where the two terms in the square brackets originate from 
the intermediate states with empty and doubly-occupied 
chain, respectively. The many-body nature of the two- 
electron transition manifests itself in the fact that these 
terms cancel each other in the limit U — ► 0. 

With the help of definitions (1) and (2), Eq. (3) can 
be presented as 



r ; ( £l - £ )r r ( £2 - £ ) 

<L 2nhe>(e + Uf ^ 



e) 1 - fr(£2 - e) 
(4) 



In principle, the integrand needs to be regularized at 
the resonances, e — 0,—U, following, e.g., the proce- 
dure of Ref. 25. However, the corresponding singular 
contributions are negligibly small in the cases of interest, 
l £ 2 S> T, where T is the temperature. 

We will establish the analytical form of t c in the 
limit of large enough bias, V, when it is temperature- 
independent. In this limit, the integration in Eq. (4) 
is restricted to the interval e\ — V/2 < e < E2 + V/2, 
which is set by the distribution functions in the leads. 
Neglecting the weak energy dependencies of Ti and T r 
and evaluating the integral, we obtain 
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The asymptotes of t c in the limits of small and large in- 
tersite interaction, U, can easily be found from 
Eqs. (5), (6) 
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It is worth noting that the rate Eq. (3) of the interaction- 
induced hop becomes [/-independent in the limit of 



strong intersite interaction, as follows from Eq. (3) as 
well as from Eq. (8). Note also that in both limits, the 
rate r c _1 is proportional to (e 2 — £1 + V), which plays the 
role of "phase volume" for two-electron transitions, as we 
will see below. 

It is seen from Eq. (3) that r c oc exp [2(D — d\ — d,2)/a], 
i.e. it is shorter than the waiting time for direct tun- 
neling, which is oc exp (2D/a). Important, however, is 
that r c can be comparable to r 2 . As we will see later, 
transport, dominated by two-electron transitions, is most 
prominent when r c and r 2 are of the same order. 



III. NON-OHMIC REGIME 
A. QUALITATIVE DISCUSSION 

In this Section we assume that the temperature is zero, 
so that hopping transport through a chain is possible if 
all hops, involved in the passage of current, are activa- 
tionless. All throughout the paper we assume the on-site 
Hubbard repulsion, Uo, to be infinite. To classify differ- 
ent transport regimes, we first note that when the bias, 
V, is big enough, V » U, the intersite repulsion, U, can 
be neglected. The current path then consists of three 
hops, I — ► 1, 1 — ► 2, and 2 — > r, in arbitrary order, in 
the sense that the only condition for a hop to occur is 
that the initial state is occupied, while the final state is 
empty. In particular, as illustrated in Fig. 1, the wait- 
ing time for the hop I — > 1 does not depend on whether 
site 2 is occupied or empty. Similarly, the waiting time 
for the hop 2 — > r is independent of the occupation of 
site 1. This, however, does not mean that the transport 
in the regime of high bias is completely uncorrelatcd. It 
only means that the correlations are of purely Hubbard 
origin, i.e. the electron in the left lead has to "wait" 
for the hop 1 — ► 2, leaving site 1 empty, to occur. This 
Hubbard-correlated regime takes place for V > V c , where 
V c is determined by the condition s± + U = V c /2. As seen 
from Fig. 1, at V = V c , site 1, shifted upward by U, due 
to site 2 being occupied, becomes aligned with the Fermi 
level in the left lead. 

As the bias is reduced below V c , the intersite repulsion 
changes radically the electron dynamics. Now the activa- 
tionless hop Z — > 1 is possible only if site 2 is empty (see 
Fig. 1). When site 2 is occupied, the transition I — > 1 
is "blocked" . Clearly, the average current drops down in 
a step- like fashion as the bias is swept through V c . It is 
also clear that below V c the current is more correlated 
than above V c . This is because the hops I — > 1, 1 — > 2, 
and 2 — > r occur in a strict succession. 

The two above regimes were discussed in the literature 
before 5,7 . However, the most nontrivial scenario of ac- 
tivationless passage of current unfolds when both transi- 
tions I — > 1 and 2 — > r are blocked. The concept of inelas- 
tic cotunneling through a dot 22 ™ 24 suggests that in this 
case the activationlcss current is due to the two-electron 
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transition (1,2) — » (1, r). However, the theory of inelas- 
tic cotunneling implies that the two-electron transition is 
immediately followed by relaxation. For a two-site chain 
this relaxation is a single-electron transition 1 — * 2. Then 
the "immediate" relaxation requires that sites 1 and 2 
are spatially close to each other (as is the case for dots). 
Note, however, that the spatial proximity of sites 1 and 
2, meaning that di, 0I2 <C D, makes the two-electron time 
t c oc exp[2(D — di — d 2 )/a] {see Eq. (3)} quite long, and 
actually comparable to the direct-tunneling time. Thus, 
in order to yield a significant current, the separation, 
(di + d 2 ), should be ~ D. On the other hand, large 
separation of sites 1 and 2 unavoidably opens up new 
relaxation channels, which are "reverse" single-electron 
transitions 1 — > I and r — ► 2. Obviously, these transi- 
tions "undo" the two-electron transition and forbid the 
current to flow. Then, it is quite nontrivial that for cer- 
tain energy configurations {51,62} of the sites, the re- 
verse relaxation channels are blocked, leaving the tran- 
sition 1 — > 2, the only allowed transition that can follow 
a two-electron hop. At the same time, the transition 
1 — ¥ 2, happening after the two-electron hop, completes 
the current cycle. We will dub such configurations with 
unidirectional relaxation as ratchet configurations. 



B. RATCHET-TYPE CONFIGURATIONS 



FIG. 2. Schematic illustration of the non-ohmic transport 
regime dominated by two-electron hops (1,2) — > (1, r) fol- 
lowed by an inelastic transition 1 — > 2. For the ratchet-type 
configuration shown, both "backwards" transitions following 
(1,2) — ► (1, r) are forbidden. 



The condition allowing a two-electron transition has 
the form Eq. (9) since this transition can be viewed as 
a transfer of an electron from the left to the right lead, 
accompanied by a "backward" excitation 2 — ► 1 of the 
other electron. Obviously, the "direct" transition 1 — * 2 
can occur after the two-electron hop. The ratchet-type 
configuration is such that this direct transition is not pro- 
ceeded by other single-electron transitions. This would 
absolutely be the case if the following two conditions are 
met 



£1 < j, (10) 
e 2 + U>~. (11) 



Consider the configuration of sites depicted in Fig. 2. 
Site 2 is lower than — V/2 and thus is occupied. The 
energy, Si + U, of site 1 is above V/2, so that site 1 is 
empty. For this configuration both transitions I — ► 1 and 
2 — > r are blocked. The condition that the two-electron 
transition (1,2) — > (l,r) is allowed at zero temperature 
reads 



The lhs in Eq. (9) is the minimal energy required to 
transfer an electron between the leads, while the rhs is 
the energy required to transfer an electron between sites 
2 and 1. 




Indeed, the first condition ensures that the electron 
from the occupied site 1 cannot hop back into the left 
lead, while the second condition guarantees that the 
empty site 2 does not get occupied as a result of an elec- 
tron hop from the right lead. Our prime observation is 
that ratchet-type configurations are possible, i.e. all five 
conditions, (ei + U) > V/2, e 2 < -V/2, and Eqs. (9)- 
(11), are satisfied simultaneously within a finite domain 
on the (si,£ 2 ) plane. This domain corresponds to the 
triangular region in Fig. 3, where I c denotes the magni- 
tude of current in the regime of two-electron hops. This 
region is situated adjacently to two rectangular domains, 
within which one of the transitions I — > 1 or 2 — > r is 
blocked. I\ stands for the magnitude of the strictly se- 
quential current in these domains. Finally, within the 
domain ( E\ < —U + V/2, e 2 > —V/2 ) both transitions 
are unblocked, so that the current, I 2 , is limited only by 
the Hubbard correlations. 
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FIG. 3. (Color online) Domains of the site energies corre- 
sponding to the different regimes of transport. 



To summarize our qualitative analysis, in Fig. 4 we 
depict schematically the current-voltage characteristics 
at zero temperature. In the absence of interactions the 
I — V curve would exhibit only one step, from the "long- 
hop" current, Id, corresponding to direct tunneling be- 
tween the leads, to Hubbard-correlated current, I 2 . In- 
teractions give rise to two additional steps: from Id to 



"doubly-blockaded" 
blockaded" current, 



current. 



I c , and from I c to "singly- 



ei-£ 2 2 led 



V. 



FIG. 4. Schematic I — V charac- 

teristics of the two-site chain at zero temperature is shown 
for the case — U < e\ + £2 < 0. Weak bias dependencies of 
current within each step are neglected. 



C. CALCULATION OF CURRENT 

For infinite Hubbard repulsion, possible sets of the 
occupation numbers (711,713) of the sites 1 and 2 are 
restricted to (0,0), (0,1), (1,0), and (1,1). Following 
Refs. 1, 4 we introduce the probabilities, P ni ,n 2 , of each 
set. Normalization requires that Po^+^o. 1+^1,0+^1.1 = 
1. 

As discussed above, in the regime V > V c , when both 
transitions I — > 1 and 2 — > r are unblocked, the trans- 
port involves all four sets of the occupation numbers. 
For V < V c , when Z — > 1 is blocked, double occupa- 
tion is prohibited, i.e. Pi 1 = 0. The master equations 
for probabilities P ni ,n 2 can be cast in a form that ac- 
counts for both cases, i.e. applies within the entire inter- 
val —V/2 < £2 < £1 < V/2, as follows 





Pi,o 


Po,o 


dt 




Tl 


rfPo.i 


Po,i 


Po,i 


dt 


T 3 




dP ,o 


Po,o 




dt 


Ti 


T 3 



p 



+ — 6(1/-^), 



&(V - V c ) + 



P 



1,0 

; 

To 



(12) 
(13) 
(14) 



where 8(x) is the step-function. It ensures that Pi t i 
drops out of the system (12)-(14) for V < V c . The ex- 
pression for current that accounts for both cases reads 



/ = 



eP i , ePii 



©(V - V c ). 



(15) 



In the stationary regime, the master equations can be 
solved for P ni ,n 2 m terms of P0.1, namely, Pqq = 
Po.m/ra, Pi, = Po,i [t 2 /t 3 + (t 2 /t0 Q(V - Vc)], Pi,i = 
Po,i { T 3/ T i) Q{V — V c ). With the normalization condition, 
Po!o + Po,i + Pi,o + Pi.i© (V - Vc) = 1, we get 



P0.1 = 
for V > V r and 



Ti + T 2 + T 3 + ^ (T 2 + T 3 ) 



Po.i = 



Ti 



T 2 + T 3 



(16) 



(17) 



for V < V c - Substituting Eqs. (16), (17) into Eq. (15) 
yields 5,7 



e (Ti + t 3 ) 



(Ti + T 3 ) 2 - T X T 3 + T 2 (Ti + T 3 ) ' 



and 



h = 



V < V c . 



V > V c , (18) 



(19) 



In the remainder of this Section we calculate the values 
I\, I2, and I c , and study the correlation characteristics 
of transport within each step. 



To find I c , we note that in the regime of current domi- 
nated by two-electron transitions only Pq^ and P^o are 
nonzero. They are related via the master equation 
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dPi,o 
dt 



Pi,o 

T 2 



Pq,i 



(20) 



and the normalization condition P ,i + Pi.o = 1- I n 
the stationary regime, these two relations yield for I c = 
Pi,o/t 2 the following expression 



L = 



(21) 



Secondly, the product t ± t 3 contains the same tunneling 
exponent as r c . These observations allow one to express 
the magnitudes of the steps as follows 



Ic 

Id 



ex P(7T 



2cosh{|[f -(d 1+ d 2 )]Y 



(22) 



The I — V characteristics of a chain contains three 
steps, see Fig. 4. Their magnitudes, I c /Id, h/Ic, an d 
Ii /Ii , depend on the spatial positions of the sites. From 
Eqs. (18)-(21) it is easy to analyze the dependence of 
these magnitudes on d\ and d 2 . First, we note that the 
product TiT 2 t 3 cx cxp(2Z?/a) does not depend on d\,d 2 . 
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(23) 
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(24) 



From Eqs. (22)-(24) it is easy to see that the magnitude 
of the first step is large, as long as di, d 2 ^> a. Concerning 
the second and the third steps, the more pronounced they 
are, the smaller the difference \d\ — d 2 \. Even for a com- 
pletely symmetric arrangement, d\ = d 2 = d, the second 
step is present only if d < D/A. Precisely at d = D/A, we 
have r 2 = r c . As d decreases, the magnitude of the sec- 
ond step grows first as exp |8(D/4 — d)/a| for D/A > 

d > D/6, and then slower, as exp for 

d < D/6. The origin of this growth is that for smaller d 
the waiting time, t c , for two-electron transitions becomes 
progressively longer than t ± , r 3 . It is also possible to con- 
clude from Eq. (24) that the magnitude of the third step, 
I 2 /h, does not exceed 4/3. The maximum magnitude 
corresponds to d\ = d 2 < D/6. 



D. CORRELATION PROPERTIES 



As mentioned in the Introduction, the advantage of 
the two-site model being exactly solvable is that it al- 
lows one to analyze the applicability of the mean-field ap- 
proach, one of the main ingredients of the theory of hop- 
ping transport. The exact solution captures correlations 
in the occupation numbers of sites, whereas within the 
mean-field description these correlations are neglected. 
For the ohmic transport, the effect of correlations on the 
average current was studied in Ref. 1. Below we examine 
the correlation properties of current for various regimes 
of non-ohmic transport. 

The average occupation numbers of the sites, (ni) and 
(n 2 ), can be expressed in terms of P ni ,n 2 as ( n i) = 



Pi,a+Pi,i and (n 2 ) = Po.i+Pi,i, respectively. Neglecting 
correlations is equivalent to setting zero the difference 

- (ni)(n 2 ) = P u - (Xo + Pi,i) (P ,i + Pi,i) 

= Po,oPi,i - Po.iPifi- (25) 

We will characterize correlations by the parameter 
k = [(nin 2 ) / (ni) (n 2 )} - 1, which, using Eq. (25), 
can be presented as 



-Pq,o-Pl,i ~ PqsPi.o 
(ni)(n 2 ) 



(26) 



so that for uncorrelated transport k = 0. It can now be 
shown that, for V > V c , the exact equations (12)-(14) 
can be cast in the following form 



d{ ni ) = _ (m)(l-(n 2 )) + r { ni )(n 2 ) 

dt T t T 2 T 2 

d{n 2 ) _ (ni)(l - (n 2 )) _ (n^) _ (ni)(n 2 ) 
dt To To r, 



(27) 



In this form, the mean-field description emerges upon ne- 
glecting the last terms in the right-hand sides. Solving 
the resulting system of equations for average occupations, 
we reproduce the result of Ref. 3 
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(28) 



From Eqs. (18) and (28) it can be seen that the ratio 
h /I 2 ¥ can be expressed as a function of a dimensionless 
parameter z — r 2 (r^ 1 + t 3 _1 ) as follows 
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F(z) = 



TMF 
1 2 



b + z + yj(b + zf -4b 
2(6 + z-l) ' 



b = 



(Ti+T 8 ) S 



T,To 



(29) 



The function F(z) is plotted in Fig. 5 for different values 
of the ratio Ti/t 3 . We see that, unlike the ohmic case 1 , 
the exact current can exceed the mean-field value if r 2 is 
sufficiently long. Moreover, for z = 1, we have I 2 — 1 2 ¥ 
for arbitrary ratio T x jT^. 




FIG. 5. Ratio of exact and mean-field currents is plot- 
ted versus a dimensionless parameter z = T2 (tj 1 + t.^ 1 ) for 
n/r3 = 0.73 (solid line); t\/t3 = 0.26 (dashed-dotted line); 
and T\/T3 = 0.13 (dashed line). 

The explanation of these facts can be obtained if we 
express the correlation parameter, k, in terms of t ± , t 2 , 
and r 3 , using Eq. (16) 



k = 



T?(l-Z) 



Ti[TiT 3 - T 2 (Ti + T 3 )] 

(ti + r 3 )(riT 2 + t 2 t 3 + r 3 2 ) (r x + t 3 )(ti2 + r 3 ) 



(30) 

We see that for z = 1 we have k = 0, so that the cur- 
rent is effectively uncorrelated. For this reason we have 
I 2 = P/ F when z = 1. For z < 1, the correlation parame- 
ter is positive, resulting in I 2 < I 2 F . Similarly, I 2 > I 2 F 
for z > 1 is the manifestation of the fact that the current 
is negatively correlated. The latter has a simple explana- 
tion: for long r 2 (i.e. for z 3> 1), as an electron tunnels 
onto site 1 from the left lead, site 2 is likely to be empty. 
Then the electron does not have to "wait" extra time be- 
yond t 2 to proceed to site 2. The mean- field description 
does not capture this effect, thus causing I 2 lF < h- 

For V < V c the mean-field description reduces to the 
following equations for the average occupation numbers 



d(ni) 
dt 



1 



<ni»[l-(n2>] 



_ ( ni )(l-(n 2 )) 

1 

T 2 

d(n 2 ) = (ni)(l - (n 2 )) _ (n 2 ) 

dt T 2 T 3 



(31) 



where the factor (1 — (n 2 )) in the first equation ex- 
presses the fact that Z — > 1 is possible only when site 
2 is empty. It easy to see that these equations yield 
7J VIF = ej (Ti + r 2 + t 3 ), i.e. the exact value of current in 
the regime V < V c , see Eq. (19). The explanation why 
neglecting the correlations does not change the average 
current can be inferred from rewriting Eqs. (31) in the 
following form 



Ti d(ni) 


r(l-(n 2 » <m>- 


! <«i>(n2> 


Tj + t 2 dt 


- Ti + T 2 


T 2 - 


T 2 


d(n 2 ) 


r («2> (ni>- 


< 


ni)(n 2 ) 


dt 


- T 3 T 2 - 




T 2 



(32) 



Note now, that if one neglects the product (n\)(n 2 ) in 
the right-hand sides, then the mean-field equations (32) 
reduce to the exact equations (12)-(14) for V < V c , since 



for V < V c we have (ni) = Pi and (n 2 ) = P^ 



o.i- 



Then 



the fact that both Eqs. (32) and (12)-(14) yield the same 
value of (n 2 ), and thus the same current I\, is a con- 
sequence of the observation that, upon adding the two 
equations (32), the products (ni)(n 2 ) cancel each other. 

As follows from Eq. (25), in two other regimes with 
I = I\ and I = I c we have k = — 1. This reflects the fact 
that the passage of current occurs via a single repeating 
cycle in both regimes. In particular, the cycle for I = I c 
consists of two-electron transition (1,2) —> (l,r) followed 
by a single-electron transition 1 — ► 2. The corresponding 
waiting times obey the Poisson distributions with aver- 
ages t c and t 2 , respectively. This allows one to find the 
Fano factor of the current noise in the regime / = I c to be 
(t 2 + t 2 2 )/(t c + t 2 ) 2 < 1. Current noise in the Hubbard- 



correlated regime, 
7. 



I 2 , was studied in detail in Rcf. 



Ti 



IV. OHMIC REGIME 

Similarly to sequential hopping, the transport regime 
dominated by two-electron transitions is activationlcss 
only if applied bias, V , is high enough. Indeed, as seen 
from Fig. 3, the domain where the current is equal to 
I c exists only for finite V. As the bias is reduced, the 
current assumes activational character in both regimes. 
Then the question arises as to whether the two-electron 
regime survives in the ohmic limit V<T, where T is the 
temperature. We will address this question in the follow- 
ing sequence. First we derive an expression, analogous 
to Eq. (21), for the ohmic resistance due to two-electron 
transitions. Then we will demonstrate that this expres- 
sion applies to the two-site chain within a certain domain 
of energies, £i,s 2 , and positions d\,d 2 of the sites. 
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In order to generalize the derivation of Eq. (21) to 
the ohmic case, we introduce, along with the time r 2 , the 
time of a "reverse" hop 



r 2 exp 



f fi - £2 \ 
\ T J" 



(33) 



We also note that in the regime V <C T the forward 
two-electron transition, illustrated with forward arrows 
in Fig. 2, acquires an activation energy, E\ — e 2 — V, 
whereas the reverse two-electron transition (1, r) — > (I, 2) 
becomes activationless with the characteristic time, r c , 
calculated as above. Therefore, the time of the forward 
two-electron transition is 



t c exp 



where r c is now given by 



£1 - £2 - V 



1 = r ; ry 



(34) 



(35) 



with $(x) defined by Eq. (6). The latter expression was 
derived assuming that (e\ — £2) 3> T. As we will see be- 
low, in the ohmic regime, the most interesting situation 
corresponds to (£1 — £2) ~ T <C U. Then one has to 
perform the integration in the general expression Eq. (4) 
[with £1 and £2 interchanged due to the redefinition of r c 
by Eq. (34)] using the explicit form of the distribution 
functions in the leads. This yields 



1 



T,r r U 2 



£1 - £2 



t c 2Khel (£1 + U) l-exp [e 2 -e l )/T 



(36) 



As in the previous section, we assume here that the 
sites are in the ratchet configuration, i.e. their energies 
are within the interval [0, U] below the Fermi level, so 
that interaction-elevated energies of both sites are above 
the Fermi level. In the regime dominated by two-electron 
transitions, the only sets of the occupation numbers in- 
volved in the passage of current are (0, 1) and (1,0). In 
the presence of reverse transitions, the master equation 
Eq. (20) assumes the form 



dPi,o 
dt 



= -Pi 




(37) 



With the normalization condition, Po,i + Pi, a = 1, this 
yields the following stationary solution for the probabil- 
ity P ,i 



^0,1 



r c - 1 +r 2 - 1 + (r c+ )- 1 + (r-)- 



(38) 



In the presence of reverse hops, the stationary current 
can be calculated, e.g., as the difference of the forward 
and reverse currents between sites 1 and 2, namely, 



h = Pi,o/t 2 - Po,i/t~. With the use of Eq. (38), we 
then obtain 



Ir = 



1 + (r 2 /r c ) + (r 2 /r+) + (t 2 /t 2 ") ' 



(39) 



In the ohmic regime, VcT, the numerator of Eq. (39) 
is proportional to V. Substituting (34) and (33) into (39), 
we find the ohmic resistance caused by two-electron hops 



V T(r 2 + r c ) 



exp 



£2 



T 



1 



(40) 



Although in our derivation we assumed that £1 > £ 2 , 
by introducing |...| into Eq. (40) and also into Eq. (36), 
we ensure that it applies to e\ < e 2 as well. Comparing 
this result to "non-ohmic" Eq. (21), we conclude that 
basically in the ohmic regime the current, I c , acquires 
the activation energy, \e\ — e 2 \- We would like to empha- 
size that the necessary condition for both Eq. (21) and 
Eq. (40) to apply is that sites 1 and 2 are in the ratchet 
configuration. 

While in the non-ohmic regime with I — > 1 and 2 — > r 
blocked, the two-electron hops constituted the only chan- 
nel for activationless current flow, it is not obvious 
whether the "two-electron resistance" Eq. (40) can deter- 
mine the net resistance in the ohmic regime. We address 
this question below. 

Note first that, whether the current is governed by se- 
quential hopping or by two-electron hops, the transfer of 
electrons between the leads necessarily includes the tran- 
sition 1 — > 2 (or 2 — > 1). Thus, the "competition" be- 
tween two-electron and sequential hopping is decided by 
the dominant mechanism of reoccupation of sites 1 and 2 
from the leads during intervals between the hops 1 — > 2 (or 
2 — > 1). If this dominant mechanism is two-electron tran- 
sitions (1,2) — ► (1, r) or reverse, then the longest waiting 
time for this transition, r c + , must be shorter than all char- 
acteristic single-electron times involved in the reoccupa- 
tion of the sites. For single-electron transitions I < — ► 1 
the shortest such time is given by 



In 7f 



di 



Ai 

T ' 



(41) 



where Ai = min {|£i|, |£i + U\} is the smallest activa- 
tion energy of the hops 1 — ► I (when site 2 is empty) and 
I — > 1 (when site 2 is occupied). Analogously, for transi- 
tions r < — ► 2, the shortest reoccupation time is defined 
as 



In 



d 2 + 



T ' 



(42) 



with A 2 = min {|e 2 |, \£2 + U\}. The times (41) and 
(42) should be compared to the longest time of the two- 
electron transition, which in terms of d\ and d 2 can be 
expressed as follows 
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(43) 



Then the conditions r+ < and r+ < rf- can be cast 
in the form 



d 2 > 
di > 



D afAi-|ei-E 2 | 



2 2T 
D a(A 2 -| £ i -e 2 
~2 



2T 



(44) 



It is seen from Eq. (44) that the most favorable sit- 
uation for two-electron transport to dominate is when 
the energies of the sites are close to each other, i.e. 
when \si — £ 2 | *C |ei|, e 2 |- In this situation we have 
Ai « A 2 = A, so that the conditions (44) can be pre- 



sented as g?i , d 2 > 



Graphically these 



D/2 - aA/2T 

conditions are represented by the vertical and horizontal 
lines on the plane {di,d 2 }, see Fig. 6. 



D/2 



D/2 - aA/2T 




D/2 - aA/2T 



D/2 



FIG. 6. (Color online) Triangular region denned by condi- 
tions di,d 2 > D/2 — aA/2T and di + d 2 < D/2 indicates the 
site positions where two-electron hops dominate the trans- 
port. 

In the domain where these conditions are met, Eq. (40) 
applies. Note, however, that when r 2 is much longer 
than t c , it is the hop 1 < — ► 2 that determines the re- 
sistance regardless of whether the reoccupation of sites 
1 and 2 occurs via two-electron or single-electron tran- 
sitions. Thus, the "true" two-electron-hopping regime is 
realized when r 2 < r c . The latter condition can be pre- 
sented as g?i + d 2 < D/2. It corresponds to the diagonal 
line in Fig. 6. This line together with the horizontal and 
vertical lines define a non-empty domain of site positions 
d\ and d 2 if the condition 



D < 



2aA 



(45) 



is met. The condition (45) has a transparent physical 
meaning. Single-electron hops require shorter tunneling 
distances, but involve activation, whereas two-electron 
hops require tunneling over longer distances, but they 
are practically activationless when e± ~ £ 2 . Thus, in 
accordance with Eq. (45), the lower the temperature, 
the wider is the range of the barrier thicknesses where 
two-electron hops dominate the transport. 

According to the definition of A as the minimal of 
the activation energies of the hop onto (or out of) 
Coulomb-shifted and unshifted sites, its value cannot ex- 
ceed U/2. Then the most favorable situation for two- 
electron-dominated hops is E\ « e 2 w —U/2. In this 
situation, Eq. (45) can be presented as D < aU/T. 



V. CONCLUDING REMARKS 

In the present paper we studied the effect of inter- 
actions on inelastic transport through a two-site chain. 
In particular, we have identified the domain of energies 
and positions of sites where the transport occurs via two- 
electron transitions, which are impossible without inter- 
actions. In the ohmic regime, the most favorable situa- 
tion for the two-electron transport regime realizes when 
the energies of the sites are close, s\ * e 2 . Note, how- 
ever, that throughout the paper we assumed that the 
difference between the site energies, as well as their en- 
ergy distance to the Fermi level, exceed the site widths, 
Ti, T r . This ensures that there is no resonant tunneling 
of a single electron via the two-site configuration 2 . While 
interaction-induced resonant tunneling with T/, T r deter- 
mined by two-electron virtual transitions 26 does not re- 
quire complete alignment of site energies, it does require 
more than two sites. 

For a thick barrier with localized states the current is 
governed by impurity chains containing many sites 27,28 . 
Their energies are close to the Fermi level, and they are 
almost equidistant. We do not expect interactions to 
have a strong effect on hopping through these "optimal" 
chains for the following reason. While the waiting times 
for all hops constituting the chain are approximately the 
same, the net resistance is determined by the longest of 
these times. This is reminiscent of the two-site chain, 
considered above, with r 2 longer than r 1; r 3 , and t c . As 
we saw before, under this condition, the fashion (corre- 
lated or non-correlated) in which electrons were delivered 
to/from the sites of the main hop was not important. 

Unlike inelastic cotunncling, the regime of elastic co- 
tunneling through a quantum dot 24 does not have an 
analog within the two-site model. This is because in the 
case of a quantum dot an electron can be transferred 
between the left and right leads only via the dot. By 
contrast, in the two-site model, an electron can tunnel 
directly between the leads; the amplitude of this process 
exceeds the amplitude of elastic cotunncling via each of 
the two sites. 







Note also, that the two-site model is too simplistic to 
capture another interaction effect on hopping transport, 
which is the effect of distant "fluctuators" away from the 
current-carrying path on the passage of current 29 . 

There exists certain analogy between the regime of 
transport by two-electron hops via impurities consid- 
ered here and transport by cotunncling through gran- 
ular arrays 30 ^ 32 . The important difference is, however, 
that in Refs. 30-32 each granule is assumed to contain 
a large number of levels; cotunneling proceeds via these 
levels, and it is the finite charging energy of an individual 
granule that allows this process. Interaction between the 
granules is essentially irrelevant in calculating the rates 
of hops 30 " 32 . 

Since we consider sites rather than granules, it is the 
on-site Hubbard interaction, Uo, that plays the role of 
charging energy. This interaction is assumed infinite in 
our consideration. On the other hand, the intersite in- 
teraction plays a central role for transport. As a con- 
sequence of the fact that two-electron transitions in a 
two-site model are due to the intersite interaction, the 
net rate of a hop is, essentially, the sum of the rates for 
two-electron and subsequent single-electron transitions. 
By contrast, in the case of granules, the probability of 
a particular hop between two distant granules is propor- 
tional to the product of tunneling probabilities between 
all adjacent granules along the path. 
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